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Yuko UrakawaQ and Takahiro Tanaka^j] 

1 Department of Physics, Waseda University, Ohkubo 3-4-1, Shinjuku, Tokyo 169-8555, Japan 
2 Yukawa Institute for Theoretical Physics, Kyoto university, Kyoto, 606-8502, Japan 

(Dated: May 24, 2011) 

It has been an issue of debate whether the inflationary infrared(IR) divergences are physical or not. Our claim 
is that, at least, in single-field models, the answer is "No," and that the spurious IR divergence is originating from 
the careless treatment of the gauge modes. In our previous work we have explicitly shown that the IR divergence 
is absent in the genuine gauge-invariant quantity at the leading order in the slow-roll approximation. We extend 
our argument to include higher-order slow-roll corrections and the contributions from the gravitational waves. 
The key issue is to assure the gauge invariance in the choice of the initial vacuum, which is a new concept that 
has not been considered in conventional calculations. 



I. INTRODUCTION 

The importance of the gauge-invariant perturbation has 
been widely recognized since decays ago, but some subtle 
issues get in the way of its realization particularly in non- 
linear perturbation theory. During inflation, massless fields 
are known to yield the scale invariant spectrum P(k) oc 1/k 3 
at linear order. These fields contribute to the one-loop diagram 
with the four point interaction as J d 3 k/k 3 in the long wave- 
length limit, which leads to the logarithmic divergence 0]- 
l20ll . In our previous work I2H1 . we pointed out that, at least 
in single field models, the IR divergences are attributed to the 
bad treatment of gauge degrees of freedom. The gauge de- 
grees of freedom can be classified into two classes: the local 
ones and the non-local ones. In the usual calculation only 
the former is fixed by adapting particular gauge conditions at 
each space-time point. However, this is not sufficient in order 
to accomplish the complete gauge fixing because of the pres- 
ence of the non-local gauge degrees of freedom, which are 
typically the degrees of freedom to specify boundary condi- 
tions in solving the lapse function and the shift vector. These 
non-local gauge degrees of freedom are formally fixed by im- 
posing the regularity at spatial infinity in the conventional per- 
turbation theory. As a result, however, the time evolution of 
so-called gauge invariant variables is affected by the informa- 
tion from infinitely large volume outside our observable re- 
gion. We claimed that this is the origin of IR divergences I2H1 . 

Along the line mentioned above, the IR divergence prob- 
lem reminds us of the importance of maintaining the gauge- 
invariance in cosmological perturbation. In our previous 
work 1I22L we provide one simple but calculable example of 
genuine gauge-invariant quantities, and showed its regularity 
at the leading order in the slow-roll approximation. In or- 
der to realize IR regular perturbation theory, one important 
additional aspect is to guarantee the gauge invariance of the 
initial quantum state. We found that by choosing the Bunch- 
Davies vacuum, which yields the scale-invariant spectrum, at 
the lowest order in slow-roll approximation, the gauge invari- 
ance of the initial quantum state is realized. In this paper, 



we extend our argument about IR regularity of such genuine 
gauge invariant quantities and the existence of gauge-invariant 
initial state to the quadratic order in the slow-roll approxima- 
tion. This extension would be wanted, because the presence 
of IR divergences that has been reported so far mostly starts 
with this order (see Ref. fl9ll for a recent review). We also in- 
clude the discussion about the contributions from the graviton 
loops. 

To quantify the primordial fluctuations and provide the 
testable predictions for models of inflation, it is necessary to 
remedy the singular behaviour of IR corrections as well as the 
ultraviolet divergence i23U25ll . The feasibility of the secular 
growth of IR contributions has also been addressed, motivated 
as a possible solution to the smallness of the cosmological 
constant lH^S. (See also Refs. JH and JH.) Despite the 
several efforts I^|32| - |43[| . the debate regarding the possibil- 
ity of the IR divergence has not been settled. To put an end 
to this debate, following the idea presented in our previous 
works I2IH22II . we explicitly show the absence of the IR di- 
vergence, restricting our argument to single field models of 
inflation. 

Our paper is organized as follows. In Sec. [Til we give the 
setup of our problem and briefly review our solution to the IR 
divergence problem. In Sec. Ill Bl we clarify the relation be- 
tween the residual gauge degrees of freedom and the boundary 
conditions in solving the lapse function and the shift vector. 
In Sec. [Ill] we give one example of genuine gauge-invariant 
variables. In Sec. [TV] we show the regularity of the gen- 
uine gauge-invariant variable and study the requirement of the 
gauge-invariance on the initial quantum state. Our results are 
summarized in Sec.lVl 
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In this section, we briefly summarize our solution to the IR 
divergence problem, proposed in our previous work l22tl . 
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A. Basic equations 

We consider a standard single field inflation model whose 
action takes the form 

M 2 r 

s = ^ I V=g[R- g^Mu - 2V(<t>)]d i x , (i) 

where M p i is the Planck mass and the scalar field <fi was 
rescaled as <fi — > 4>/M p \ to be dimensionless. The ADM for- 
malism has been utilized to derive the action of the dynam- 
ical variables particularly in the non-linear perturbation the- 
ory l44ll . Using the decomposed metric 

ds 2 = -N 2 dt 2 + hijidx 1 + N l dt)(dx*> + N j dt) , (2) 

the action is rewritten as 



S 



N S R~ 2NV(4>) + —{E irj E ij - E 2 ) 



N l d l( /)) 2 - XI," O.oi),*:* 



dV 



(3) 



where S R is the three-dimensional scalar curvature and E, 
and E are defined by 

1 



Ei 



(d t hi 



E = h v Ei 



(4) 



The spatial index of Ni is raised and lowered by hij . 

In this paper we work both in the comoving gauge and in 
the flat gauge. We defer the introduction of the latter gauge to 
Sec. IIH Bl The comoving gauge is defined by 



= 



(5) 



where 8<fi is the perturbation of the scalar field. We decompose 
the spatial metric as 



(6) 



where a := e p denotes the background scale factor, and 
tr[5-f] = 0. Using the degrees of freedom in the choice of 
the spatial coordinates, we further impose the gauge condi- 
tions l 8jij = 0. Here, the indices of spatial derivatives are 
raised or lowered by using Kronecker's delta as l = S l; >0j. 

Varying the action with respect to N and N l , we obtain the 
Hamiltonian and momentum constraints as 



S R-2V - N- 2 {E l3 E lJ - E 2 ) - N~ 2 [d t 



= 0, (7) 
(8) 



Dj [AT 1 =0. 
Introducing the perturbed variables as 

:= e- 2p h l3 , N, = e p N t , N l := h ij Ni = e p N l , (9) 
we factorize the scale factor from the metric as 



ds 2 



(N 2 - N.N^dr] 2 

+ 2N t dr]dx l + h,,d.r'd.r' . (10) 



Expanding the perturbations, Q = SN(:= N — 1), Ni, £, and 
S^ij as Q = Qi + Q2 + • • ■ , the zeroth-order Hamiltonian con- 
straint equation yields the background Friedmann equation: 



6p' 2 = <p' 2 + 2e 2p V(4>) , 



(11) 



where a prime " ' " denotes the differentiation with respect to 
the conformal time 77. The constraint equations at the linear 
order are obtained as 

e^VSNt - + d 2 Ci + i>'0'.\,.; = , (12) 
49, (p'SNt - CO - 2 N til + 0,0*. \ r : = , (13) 

where d 2 := l 0i. The higher-order constraints can be ob- 
tained similarly. 



B. Residual gauge degrees of freedom 

The constraint equations (O and dH) allow us to describe 
the non-dynamical variables N and Ni in terms of £. Here 
we stress that the constraints ( fT2b and the divergence of ( fT3] l 
are elliptic-type equations, which require boundary conditions 
to solve. Even though we impose the gauge conditions (0 
and (O at each space-time point, N and Ni are not uniquely 
determined because of the presence of such non-local gauge 
degrees of freedom. At the first order of perturbation, these 
degrees of freedom are studied in Ref. 12211 . where general 
solutions of SNi and iV^i are given in the form: 



1 



(14) 



Here, an arbitrary vector function Gi(x) that satisfies the 
Laplace equation d 2 Gi(x) = was introduced to make ex- 
plicit the presence of degrees of freedom corresponding to the 
boundary conditions. Substituting Eqs. ( fT4b and (TT~5T > into the 
equations of motion for fi and £7^,1, we find that the intro- 
duction of the gauge function Gi(x) modifies their evolution 
equations as well 112211 . 

The ambiguity originating from the choice of the vector 
Gi(x) is a sign of the presence of residual gauge degrees of 
freedom. Here, we explicitly show that Gi(x) represents the 
residual gauge degrees of freedom that remain undetermined 
even after specifying the gauge by the conditions © and (O. 
Since the gauge condition 8(f) = completely fixes the tem- 
poral gauge, the residual gauge can reside only in changing 
the spatial coordinates: x % — > x l = x l + 8x % . The metric 
perturbations then transform as 

#i,iO) = Ni, x {x) - 5x'i , (16) 
Ci(x) = Ci(:c) - \&&Xi , (17) 
Sjij,i(x) = 6jij t i(x) - 2 (d(i5xj) - \d k 5xk8ij) . (18) 
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In this section we associate a tilde " " " with the perturbed 
variables in the gauge with Gi ^ to discriminate them from 
the perturbed variables in the gauge with G; = 0. 

Since we have not changed the temporal coordinate, the 
lapse function remains unchanged at the linear order. Equat- 
ing 5N\ with SNi, given by Eq. (TBI , we find that (i is related 
to (i as 



Ci = Ci + \d l G % . 
Comparing Eq. (T7\ with (T% , we obtain 



(19) 



(20) 



Imposing the transverse condition on Sjij in Eq. ( fTST l. we ob- 
tain another condition for 5x l as 

d 2 8 Xl = -Idi&Sxj . (21) 

Then, Eqs. (l20l i and (fSTJ are integrated to give 



dr)Gi(x) + - / drfdid^&iGjix) 



+ dr)hi(x) + Hi(x), 



(22) 



where we introduced vector functions hi(x) and Hi(x) that 
satisfy 



(23) 
(24) 



d*h l {x) = d 2 h l {x) = Q, 
3d 2 H i (x)+d i d j H j (x) = . 

Using Eqs. ( TT7b and ( l20b . Eq. (TT~5T > is recast into 

tii,i{x) = N itl (x) + ^di&Sxj 
5p 

_ ldid- 2 d j G jA (x) + G iA (x) . (25) 

Comparing Eq. d25l l with Eq. ( fT6] l, the vector function hi(x) 
is determined by 

h i( x ) = T7 / dridi^Gjix) + ^H^x) . (26) 

The degrees of freedom in boundary conditions are then found 
to represent the change of the spatial coordinates: 



Sxi(x) = - jdT 1 G l (x) + ^ ldr]did- 2 d J G 3 {x) 



1 

4 

+ \Jj J dj l di&G^x) + Hi(x) 



(27) 



which is basically expressed in terms of Gi(x). We note that 
the time-independent vector Hi{x) can be absorbed into the 
integration constant of the temporal integral, J dr/ {— Gi(x) + 
(l/4)a 4 a- 2 a j G J (a;)}. Substituting Eq. <E7) into Eqs. CE3 



and (TT~8T >. we find that spatial components of metric perturba- 
tion transform as 



Ci 



C1 + l/' 



1 



a 4 a 



(28) 



5%j,i = + dri {2d {l G n - \(d t d,jd~ 2 + 5 l0 )d k G k ) 



I jd 2 d (l H j} . 



(29) 



When we consider the universe with infinite volume and re- 
quire that all quantities are regular at the spatial infinity, the 
solution of 5Ni and iVj i would be specified uniquely. How- 
ever, in this case, we observe the singular behaviour in the 
loo p co rrections of the curvature perturbation ( iH-HH [13- 
I20I l42ll . This is because the IR fluctuation acausally propa- 
gates through the non-physical gauge modes and comes into 
play. In contrast, if we do not care about any singular behav- 
iors at infinity, that would never be observed by us, a variety 
of homogeneous solutions Gi can be added to the solution of 
8N\ and JVi 1. In our previous work 12111 . we have shown that, 
by choosing the function Gi appropriately, we can guarantee 
the regularity of fluctuations of Q in the flat gauge as long as 
a finite spatial region of our universe is concerned. We think 
that this is a remarkable progress, but the prescription given 
in Ref. I2H1 is not completely satisfactory in that the loop cor- 
rections for £ depend on the choice of the boundary conditions 
for the lapse function and the shift vector. This fact signifies 
that the n-point functions for the "so-called" curvature per- 
turbation ( is not a genuine gauge-invariant quantity, when 
we take into account the gauge degrees of freedom associated 
with the choice of boundary conditions. Our discussion here 
can be extended straightforwardly to the higher-order in per- 
turbations. 



m. GAUGE-INVARIANT QUANTITIES 

In this section, we provide one simple example of genuine 
gauge-invariant quantities. If we compute genuine gauge- 
invariant quantities, the results by definition should be unaf- 
fected by the choice of the gauge. Hence, they should be IR 
regular even if we calculate them based on the standard per- 
turbation theory. We demonstrate this in the following two 
sections. 



A. Definitions of scalar curvatures 

One simple way to realize the gauge invariance is to use 
variables defined in a completely fixed slicing and threading. 
What is revealed in the previous section is the fact that the gen- 
uine gauge-invariant variables cannot be constructed by sim- 
ply adapting gauge conditions to metric components at each 
space-time point. In order to fix the gauge completely, we 
also need to remove the unphysical degrees of freedom asso- 
ciated with the choice of boundary conditions. This cannot 
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be achieved easily due to the difficulties in removing all arbi- 
trariness regarding the choice of space-time coordinates. It is, 
however, possible to calculate genuine gauge-invariant quan- 
tities even if we do not accomplish the complete gauge fixing. 

Since the time slicing is uniquely fixed by the gauge con- 
dition 6<f> = 0, it is enough if we can arrange quantities so 
as to be invariant under the transformation of spatial coordi- 
nates. In our previous work 12211 . we proposed to calculate 
n-point functions for the scalar curvature of the induced met- 
ric on a <f) =constant surface, S R. Although S R itself does not 
remain invariant but transforms as a scalar quantity under the 
change of spatial coordinates, the gauge invariance of the n- 
point functions of S R would be ensured, if we could specify 
its n arguments in a coordinate-independent manner. The dis- 
tances of spatial geodesies that connect pairs of n points char- 
acterize the configuration in a coordinate independent manner. 
Based on this idea, we specify the n spatial points in terms of 
the geodesic distances and the directional cosines, measured 
from a reference point. Although we cannot specify the ref- 
erence point and frame in a coordinate independent manner, 
this gauge dependence would not matter as long as we are in- 
terested in the correlation functions in a quantum state that 
respects the spatial homogeneity and isotropy of the universe. 

We consider the three-dimensional geodesies whose affine 
parameter ranges from A = to 1 with the initial "velocity" 
given by 



(1A 



number of the contained interaction picture field operators. 
For simplicity, we neglect the terms that do not yield the IR 
divergence. In the above expression, each term contains two 
pairs of contraction between creation and annihilation opera- 
tors. Only when one of these pairs does not contain any dif- 
ferentiation, the term potentially contributes to IR divergence. 
Since the loop integrals diverge at most logarithmically, one 
spatial or temporal derivative is sufficient to remedy their di- 
vergent behaviors. Noticing that the curvature perturbation 
in the first-order scalar curvature is multiplied by the spatial 
derivatives as 



(33) 



X 1 



A=0 



the terms in 9 Rs that include more than one interaction picture 
field operators with spatial or temporal derivatives do not yield 
divergences. This statement also applies to the terms in 9 i?2. 
Since 9 i?2 contains at least one interaction picture field oper- 
ator that is differentiated, the terms in 9 i?2 that include more 
than one differentiated interaction picture field operators do 
not yield IR divergences. 

The loops of gravitational wave perturbation without 
derivatives yield the logarithmic divergence, too. However, 
the gravitational wave perturbation 6%j with derivatives no 
longer contributes to such divergent loop corrections. 

Hereafter, we denote an equality which is valid only when 

IR 

we neglect the terms irrelevant to IR divergences by « . 
Then, abbreviating the unimportant pre-factor, we simply de- 
note the scalar curvature S R as 



We identify a point in the geodesic normal coordinates X % 
with the end point of the geodesic x l (X, A = 1). Noticing 
that in the absence of the fluctuations X 1 coincides with x 1 , 
we expand x l (X) as 

x*(X) :=X t + 8x\X) . 

We denote the spatial curvature whose argument is specified 
by the geodesic normal coordinates X 1 as 



°R(ri, X) := S R( V , x\X)) 
Then, 9 R can be expanded as 



(30) 



9 R(v, X) 



Sx" 11 ■ ■ ■ Sa 



n=0 



■d in s R(r,, x% i= 



(31) 



The n-point functions of g R would be surely gauge-invariant, 
unless the initial quantum state breaks the gauge-invariance. 

Our main purpose of this paper is to demonstrate the ab- 
sence of IR divergence in the genuine gauge-invariant quanti- 
ties at one-loop order. At this order, the following three terms 
contribute to the two-point function: 



( 9 R 9 R) 4 := ( 9 i?i 9 i? 3 ) + ( 9 i?2 9 i?2> + ( 9 R 3 9 Ri) 



(32) 



where the subscripts 1, 2, 3, 4 mean the numbers of the con- 
tained creation and annihilation operators or equivalently the 



2 < [e-^fd^C 



(34) 



B. Gauge transformation 

To calculate the non-linear corrections under the slow-roll 
approximation, it is convenient to temporally work in the flat 
gauge: 

hij = e 2 " [e%. , tr[<57] = 0, 9^=0, (35) 

because all the interaction vertexes are explicitly suppressed 
by the slow -roll parameters in this gauge l2ll l44Tl . Here in 
this section we associate a tilde with the metric perturbations 
in the flat gauge to discriminate those in the comoving gauge. 
The action in this gauge is given by 



S 



IR M pi I 



TV" 1 
2Ne 
AT 



2p J2 



' + ip'- N l d z <p 



tp m -Nh^ditpdjtp 



4p'd x N'' 



dr]d 3 x. (36) 



The transformation formulae between the comoving gauge 
and the flat gauge are studied in Ref. ll-iill . and we briefly sum- 
marize them in Appendix|A] The curvature perturbation in the 
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comoving gauge £ is related to the fluctuation of the dimen- 
sionless scalar field (divided by M p i) in the flat gauge <p as 



Substituting Eqs. d42t and ( |43T > into Eq. d39l l, the evolution 
equation of £„ is recast into a rather compact expression, 



+ 3£2C f pCn +^e 2 (e 2 + 2e 3 )Cl (37) 

where we have introduced £„ := — (p'/<j)')ip, following 
Ref. 14411 . We use the horizon flow function: 



Jo 



Hi 
H 



1 d£ 



--m+1 



dp 



for to > 0, (38) 



where H is the Hubble parameter and Hi is the one at the ini- 
tial time. The horizon flow function is related to the conven- 
tional slow -roll parameters as shown in Ref. j45tl . Hereafter, 
assuming that the horizon flow functions e m with m > 1 are 
all small of 0(e), we neglect the terms of 0(e 3 ). In Eq. d37l i. 
we neglected the cubic terms that include only one graviton 
field Sjij, for the following reason. Since 9 R\ includes only 
(i, the terms in 9 R 3 that include only one graviton field S-jij, 
does not contribute to ( 9 R a R)4 after taking the contraction. 

In line with the preceding papers ll2lll32ll . in order to cal- 
culate the n-point functions, we solve the evolution equation 
(Heisenberg equation) for the operator ip, and we express ip in 
terms of the interaction picture field ipj. Variation of the total 
action with respect to ip yields 



P 2P 



Ey{m+l) 
— -i— V" 
m! 

m=0 



IR 



b' + (ff) jd t N* - N [ e "^] y didjy 5 



(39) 



where lA m ) := d m V/d(j) m . To address the regularity of the 
graviton loops, we also include the contributions from the 
gravitational wave perturbation. Variations with respect to 
the lapse function and the shift vector, respectively, yield the 
Hamiltonian constraint: 



(iV 2 - \)e 2p V + N 2 e 2p ^ -ip 



+ 2p'd^ + cb'p' + \p' 2 « , (40) 



and the momentum constraints: 



yp 

2p'd l N - N((b'd l p + d l pp) « . 



(41) 



For the calculation of one loop corrections, it is enough to 
solve the constraint equations up to the quadratic order. These 
constraint equations are solved to give 



iR e x 



difr « exC - lexe 2 CnC 



(42) 
(43) 



IB 



-2exC n + ^e 1 (4e 1 +e 2 )C 



P 



where the differential operator £ is defined by 



L := d 2 p + (3 - £l 



(44) 



(45) 



We expand ( n as ( n = C„,i + Cn,2 + Cn,3 H and denote 

£n,i simply as ip := £n,i- The equation of motion d44l is 
expanded as 



£^ = 
in 



(46) 



(47) 



IR 



n . 3 



£l£ 2 1pd p 1p - -£ 2 £ 3 i> 2 

- 2e 1 i>\ 1 d 2 i) - -fiSfa didjif) , 
p' 2 p' 2 

42 £ l(^ 2 C«,2+Cn,2dV) 



2p' p' 



2p' 



(48) 



Here, we neglected on the right-hand side of Eq. d48| ), 

because a particular solution of 57^,2 is associated with 
derivatives, and we set its homogeneous solution to zero. At 
the second order, Eq. d47b is integrated to give 



« Cn,2 + -Sf^Xidjlp, 



(49) 



where, for a later use, we have distinguished the part contain- 
ing the contributions due to gravitational waves from the pure 
scalar part £„ :2 given by 

Cn,2 « (y + £2) V- 2 + eiV^V + £i(£i + e 2 )ipd p il) 

+ 5( n , 2 + \ 2 ij(d p ~x t d l )iP . (50) 
Here, („ 2 includes the non-local term: 

"3 



-e 2 {2ex+£ 3 W 



(51) 



It should be emphasized that the homogeneous solutions £ 2 i 1 " 1 
and \ 2 ip(d p — x l di)ip can be added to C, n ,2, where the time 
dependent functions £ 2 and A 2 should be of 0(e 2 ) and their 
derivatives should be of 0(e 3 ). One can easily check that the 
above solution satisfies Eq. d47l i to the present order of ap- 
proximation, using the commutation relations 

[L,d p ] = -2(l-ex)^d 2 + 0(e 2 ), 
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[£>, x%] = -yd 2 , (52) 
Z, II A =-^(2d p + l) + 0(e). 



We are also allowed to change the solution of C, n .2 at 0(e) 
by adjusting its solution at 0(e 2 ) . In the succeeding section, 
we will explain that the solution of ^ is restricted by the 
requirement that the canonical commutation relation should 
be consistently satisfied. This requirement is, however, not 
enough to determine ( n ,2 uniquely. Therefore, in Eq. ( l50b . 
we fixed the terms of 0(e), requesting that, in addition to the 
consistency of the commutation relation, Q ni % should be kept 
in the simplest form. 

At the third-order of perturbation, Eq. d48l l is integrated to 
give 

IR " 1 1 

Cn,3 ~ Cn,3 + g (^Tl^Tl)^ Xidjtp + -5j[ 3 Sj* 1 XjXldid k lp, 

(53) 

where C„ j3 is the part purely composed of the scalar perturba- 
tion as 

Cn,3 « + A 3 V 2 (^ - X % di)ll> 



1 



1 



+ \eWd 2 ^ + -ei(6ei - e 2 )ip 2 x i d i ip. (54) 

It is again allowed to add homogeneous solutions whose coef- 
ficients £3 and A3 are of 0(e 2 ) and their derivatives by p is of 
0(e 3 ). Here we note that [L, d p - x l d l ] = 0(e). 



C. Consistency of commutation relations 

We have to take into account the following additional con- 
ditions that determine the choice of the homogeneous solution 
in £„. Until Eq. d82l . we will leave £3 unspecified, but the 
other time dependent functions £2, A 2 and A3 are constrained 
in principle so as to guarantee the normal commutation rela- 
tion for ip, as we will explain soon below. 

The evolution of the Heisenberg field £ n is usually sup- 
posed to be solved with the initial conditions that the Heisen- 
berg field £„ is identified with the interaction picture field ip 
at the initial time. This procedure guarantees that the operator 
U that relates ip to £„ by ip = U ( n U^ is unitary. In this case, 
the canonical commutation relation for ( n 



[(n(v, x ), ^n(v,v)} = iS {3) (x - y), 



(55) 



is equivalent to the commutation relation for the interaction 
picture fields 



[ip(r),x), 7r v ,(r/,y)] = iS (3) (a 



y 



(56) 



where 7r„ is the conjugate momentum of £„ and tt^, is its linear 
truncation. We give a more explicit expression only up to (e) 
relative to the leading term here. In this approximation, using 
Eq. (O, we obtain the kinetic term in the action as 



Skin = / d/y 



M, 2 ,e 2p 



dV" pl 



2N 



6' + ^') 2 + 



1 Actually, we could add a term of 0(e) proportional to i/i(dp — x 1 d{)tp to 
Cn,2. 



d J d V J d 3 xMyp £l c 

x [(1 + eiCOC + f/(e 2 - 2e 1 )C n + 0(e 2 )] . 
From this expression, we can define the conjugate momentum 

(5 Skin 

= MpVxe 2 " [2(1 + e.CnK + f/& ~ ^i)Cn + 0(e 2 )] . 

In the preceding subsection, we gave the non-linear solution 
by integrating the equation of motion without care about its 
initial conditions. Therefore unitary relation between ip and 
( n is not guaranteed. Once we obtain the definite expansions 
of and 7r„ in terms of ip and ir^,, it would be possible to 
check whether the commutation relation of ip and ir^, is guar- 
anteed from that of £ n and ir n or vise versa. We can here check 
this consistency of these commutation relations to the only 
limited extent because we have neglected the terms containing 
more than two interaction picture field operators with space- 
time differentiation. Under this limitation, we can evaluate the 
commutator, assuming [ip(rj, x), ^(77, y)} = iS^ 3 \x — y), as 

[Cn(V, x ), Kn(V,y)] 

= iS^(x -y)[l+ {46 - 3A 2 + 0(e 2 )} C„ 

+ {66 + 3A 3 + 0(e 2 )}C 2 l + •••], (57) 

where the ellipsis represents the terms containing 7Tw, and spa- 
tial derivatives, which are the beyond the scope of the present 
paper. In this way we can verify that the solution we gave is 
consistent with the expected commutation relation at 0(e). If 
we have chosen an inappropriate solution for ( n ,2 at 0(e), the 
commutation relation would not be satisfied. 

At 0(e 2 ) the terms that we could evaluate in 07] ) include 
the unspecified functions £ 2 > ^2, an d A3. From the require- 
ment that the right-hand side of Eq. d57| i should be equated to 
Eq. ( 1531 ). A2 and A 3 are related to £ 2 and £ 3 > respectively. 



D. Calculations of scalar curvature 



In this subsection, we give the expansion of the scalar cur- 
vature of a =constant hypersurface 9 R in terms of ip, which 
is the interaction picture field of £„. Using Eqs. d49l i and d53l ), 
we first perturb the scalar curvature S R, given by Eq. ( f34t . to 



7 



obtain 



s r 2 « s i? 2 + ^a^-aV, 



where and S R% are defined by 



1 



(58) 
(59) 

(60) 



It would be appropriate to emphasize that, in contrast to the 
contributions from £, the contribution from the gravitational 
wave perturbation S^ij completely cancels with each other in 
9 R. This clearly shows that the graviton loop does not lead 
to IR divergence at the one-loop order. This is essentially be- 
cause the effect on S R from IR modes of gravitational wave 
perturbation is simply caused by the associated deformation 
of the spatial coordinates. Such gauge artifacts should com- 
pletely disappear when we consider coordinate independent 
quantities like 9 R. 



s R-2 « d 2 Cn,2 + 1> d 2 (d p - 2 + e 2 /2) V, (61) 
S R 3 « 5 2 C„, 3 + L,2d 2 (d p -2 + ea/2) d 2 ^ 

+ i>{d p -2 + ea /2) d 2 ( n , 2 + d 2 ^d p ( n , 2 

-)2 



r d 2 - 4(d p 1) + |e 2 9 p 



3^2 + ^£2(^2 + £3) 



i/>. (62) 



To derive Eqs. i ]59l ) and d601> . we have used the fact that, as 
presented in Eq. ( IA8I ). the gravitational wave perturbation in 
the comoving gauge 5-jij is identical to that in the flat gauge 
Sjij, besides the terms irrelevant to IR divergences. 

Noticing that the spatial metric after removing the common 
scale factor is given by 



dA 2 = e 2C [e Sl ]..Ax l dx 3 



the geodesic normal coordinates X 1 is given by 

IR 



(63) 



(64) 



where we again abbreviated the terms that include space-time 
derivatives. The difference between the global coordinates 
and the geodesic normal ones is given by 



where 



8x l := x l (X) - X i = 5x\ + bx\ 



TD 1 

5x\ « - -StfXi , 

TD 11 

54«-( 2 X i + - 1 /, 2 X i + -(5^)»X j 



(65) 

(66) 
(67) 



Now, we are ready to calculate the scalar curvature 9 R. Sub- 
stituting Eqs. d59b and d66l l into Eq. OIK we obtain 



IR 



9 R 2 = S R 2 + SxldfR-L « S R 2 - ipX^SRi , (68) 
and substituting Eqs. d60l ) and d67l i into Eq. OIK we obtain 

1 



9 Rs 



S R 3 + 6x\di s R 2 + 8x\di s Ri + ^felfe^d/i?! 
« s i? 3 - ^ft-fia - (Cn,a + jV^ 2 ) 

+^ 2 (^^) 2 ^i (69) 



IV. IR REGULARITY AND GAUGE-INVARIANT VACUUM 

This section is devoted to show how the possibly divergent 
terms are cancelled in the n-point functions of 9 R. As de- 
scribed in Sec. IIII Al in the standard cosmological perturba- 
tion the Hilbert space has not been reduced to the one that is 
composed only of the physical degrees of freedom. Namely, 
a part of gauge degrees of freedom are left unfixed. Hence, 
an arbitrary quantum state defined in this Hilbert space can 
be non-invariant along the gauge orbit of these residual gauge 
degrees of freedom. To ensure the gauge-invariance of the tl- 
point functions of 9 R, it turns out to be crucial to set the initial 
quantum state to be gauge invariant as well. Otherwise, the n- 
point functions fail to be regular due to the gauge artifacts. 
In this section, we reveal how the gauge-invariance condi- 
tion(=regularity condition for the n-point functions) restricts 
the initial vacuum, particularly considering the one-loop cor- 
rection to the two-point function. 



A. Proof of IR regularity 

In the previous section, we expanded the scalar curvature 
9 R in terms of the interaction picture field ip. We expand ip as 



(70) 



where the creation and annihilation operators satisfy 
[ofc, cb\.i] = S^(k — k'). Focusing on the contribution from 
each Fourier mode t/j k = ffee* , the derivative operator x l di 
is rewritten as 

x l di^ k = v k d logk e lkx = d logk tp - e ikx d Xogk v k . (71) 

For illustrative purpose, we first consider the leading order 
in the slow -roll approximation j22tl . For the Bunch-Davies 
vacuum, the mode function v k is given by 

^)=-^^j- 2 ^e-^(l + ikrj), (12) 

and is easily checked to satisfy 

{d p - x%)^ k = -D k it> k , (73) 
where the operator D k is defined by 

D k :=d logk + l. (74) 
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Using Eq. ( |73| ). g R 2 and 9 R 3 could be compactly written as 

ir 1 



"R, 



IR 



-4>d D k ?p k 



9 R, 



-y&Dtrtu . (75) 



Taking the contractions of 9 R, we obtain 

(»Jk(Xi)»Ri(X 2 )} 



> 2 > 



n 



d 3 fc. t 
(2tt) 3 / 2 

x fc 2 V? (X 2 )5( 3 ) (fei-fca) 



d(logfc) 
2tt 2 



and 



IR 



n 



i=1.2 



d 3 fc, 



(27r) 3 / 2 



x (D k2 k 2 2 r k2 {X 2 )) 6^ (fci-fc a ) 

= (^)/^#5lo gfe {fc 7/2 V^(^l)}9 logfc {fc 7 / 2 ^(X 2 )}, 

where using the geodesic normal coordinate we defined 
I m := (77, X m ) for m = 1,2. Gathering the three terms 
on the right-hand side of Eq. d32l . the two-point function at 
one-loop order is summarized as 



({ 9 R(X 1 ), 9 R(X 2 )}) 4 



IR 1 



d(logfc) 



2tt 2 



5 1 2 ogfc {fc 7 ^(X 1 )^(X 2 )} + (c 



(76) 



where we symmetrized about X\ and X 2 , This indicates that 
all the potentially divergent pieces become the total derivative 
with respect to k and hence they vanish 2 . 

At the leading order in the slow-roll approximation, the 
condition d73l ). satisfied in the scale-invariant/Bunch Davies 
vacuum, was crucial to remove the IR divergences. If we do 
not choose this vacuum, the quantum state is not invariant un- 
der the residual gauge transformation, and hence the two-point 
function diverges. We think that this possible divergence is at- 
tributed to infinitely large fluctuation in the residual gauge de- 
gree of freedom corresponding to the overall rescaling of the 
spatial coordinates. This unphysical degree of freedom can 
be tamed if and probably only if we set the initial state to be 
invariant under this residual gauge transformation, as we have 
anticipated earlier. 

Now, we extend our argument to 0(e 2 ). Once we include 
the slow-roll corrections, the condition ( l73l l no longer ensures 



Here, we assumed that the ultraviolet contribution has already been regu- 
larized appropriately, say, by dimensional regularization. (See Ref. fU.) 



the gauge invariance of the initial state. Using Eqs. doTT l. d66t . 
and (|68t , the second-order scalar curvature is summarized as 



m 2 « va 2 



(1 + ei + e\ + + A 2 ) dpip - (1 + \ 2 )x i d i tp 
+ Y + 2&) V'-^- 1 £2(2ei+£3> 



Here we used 

)2 



TO Q 

&'5C„, 2 « -- [^d 2 L- 1 e 2 (2e 1 + e a )i>] , (78) 



which follows from 

£d 2 6{ n>2 « -|c [V^ 2 £- 1 e 2 (2e 1 + e 3 )^] . 

As a natural extension of the condition d73l i, we assume that 
there is a set of mode functions which satisfies 

(1 + ei + e\ + eiea) d p v k +( £l + £ -l + 2£ 2 ) v k 



-£ / 7 1 e 2 (2e 1 + e 3 )v k = -D k v k + 0(e 3 ) , (79) 



where we replaced L with L k : 



L k := d 2 p + (3 - s x + e 2 )d p + -^k 2 

P 



(80) 



We will show the presence of such mode function in the suc- 
ceeding subsection. 

For such mode functions, the second-order scalar curvature 
9 R 2 is then simply rewritten as 

-(i + \ 2 )4>d 2 D k ij 

as in the leading-order of the slow-roll approximation. Using 
Eqs. (|62]i, (|67j, and (|69]l together with Eq. (|79j, the straight- 
forward but lengthy calculation leads to 9 R 3 in a rather simple 
expression: 



9R 3 S I^ 2 



(1 + 2X 2 )D 2 k ^ - fiD k ^ 

+ (-2e\ + |eie a + \e 2 e 3 + 6£ 3 ) ip 



8Q n , 2 d 2 D k 4, 



(81) 



where we defined \i := £i + |£2-3£ 2 + i£i£2+2(£ 2 +A 3 ). To 
ensure the absence of the IR divergences, the arbitrary time- 
dependent function £3 should be chosen as 



& := \z\ - ^&2(3ei +e 3 ) , 



(82) 



to find 



"Ra « \^ 2 d 2 [(1 + 2X 2 )D 2 i; - fiD^] - 5( n , 2 d 2 D k ^ 



(83) 
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The possibly divergent terms are then summarized as 

({mm, 9 R{x 2 )}) i 

3 \m J + 2A 2 )9f ogfe - ^ ogfe } 

x{(fcVfe(X 1 )^(X 2 ))+(c.c.)} 

- (SU,2) J ^^d iogk {(k 7 MXlW k (X2))+ (CC.)} , 

(84) 

indicating that they are completely cancelled. The conditions 
on the initial vacuum state are derived by requesting the reg- 
ularity of the IR corrections. Since the IR divergence is, in 
single field models of inflation, originating from the residual 
gauge degrees of freedom l2ll I22I1 . the regularity conditions 
can be considered as the necessary condition for the gauge 
invariance. 



B. Gauge-invariant initial vacuum 

The requirement of the gauge invariance in the initial vac- 
uum leads to the condition (|79l on the mode function v k ■ Tak- 
ing the mode function v k in a similar form to Eq. ( 172b as 

«*(P) = ^^A(P). (85) 
the condition d79l can be recast into a rather simple form: 

(dp + d logk ) } k (p) + (2^ -e\- \jie2^j fk(p) 

~Lfe 2 (2e 1 + s 3 )fk(p) = 0, (86) 
where we changed the time variable p into 

p = lo g(0 '- £l + 0(e 2 ). (87) 

Notice that p is approximately identical to p in the sense 
dp/dp = 1 + 0(e). 

The mode equation L k v k = yields the evolution equation 
of f k as 

lkh(p) = Q, (88) 

where 

l k = a| + 3a p - + e- 2 ^- logfe ) 

-3(ei+e 2 /2) + 0(e 2 ). (89) 

This operator L k is identical to Lk at the leading order in the 
slow roll approximation. To fix the form of the terms written 
as 0(e 2 ) in the above equation, we need to specify the form of 
p up to 0(e 2 ). Since the explicit forms of the terms of 0(e 2 ) 
in Eq. d89l are not necessary for the following discussion, we 
leave the higher order corrections to p unspecified here. Oper- 
ating dp + 5io g A; on Eq. (|88] l, the solution of the mode equation 
( |88l is found to satisfy 

£ fc (dp + d logk ) f k = ^e 2 (2£i + e 3 )f k , (90) 



where we used [dp + d\ ogk , £fc] = -3e 2 (ei + £3/2). Mul- 
tiplying the inverse of L k on Eq. d90b , Eq. ( f90b reproduces 
the gauge-invariance condition (l86l l. where the second term 
of Eq. d86b appears as a homogeneous solution of l k . This 
indicates that the solution of the mode equation d88l consis- 
tently satisfies the gauge-invariance condition. 

It is also possible to show that the gauge-invariance con- 
dition ( f86l > is sufficient to ensure that the mode equation is 
satisfied for all wavenumbers, if it is satisfied for a particular 
wavenumber fco : &k fk (p) — 0. In fact, by using the gauge- 
invariance condition, one can show 

dl o gk£kfk\k=k = (dp + di og k) &kfk\k=k 

= lk (dp + diogfe) fk\k=k - ^ (2£l + e3) /fco + 0(e3) 

= 0(e 3 ) , (91) 

which proves that thanks to the gauge-invariance condition, 
the mode function for the another wavenumber is guaranteed 
from that for fco. 

As was anticipated in Sec. lIV Al the commutation relation 
for tjj is now verified. This commutation relation is equivalent 
to the normalization condition for v k given by 

j. ._ ie 3il (fkdpfZ - f£d p h) _ 
k ' A: 3 {l-2ei + 0(e 2 )} 

As in the case of mode equation, we assume that the normal- 
ization condition is satisfied for a particular wavelength ko as 
Nfc = 1. Then, using the gauge-invariance condition d86l ). 
we obtain 

<9io g k^k I k=k = (dp + d\ og k ) Nk I k=k 

= B-^ 2 . (92) 

where the first term B is of 0(e 2 ) and its explicit form is not 
necessary in the current discussion. An important fact is 

dpd Xogk J4 k \ k=ka =0, (93) 

holds exactly, because dp and d\ ogk commute with each other 
and the normalization condition N k is conserved. Therefore 
the right hand side of Eq. d92l ) is guaranteed to be constant in 
time. Then, by choosing £ 2 appropriately, we can always set 

d log k^k\k=k = 0(e 3 ) . (94) 

This proves that one can extend the mode function by the 
gauge-invariance condition d86l l to the other wavenumbers 
keeping the normalization condition satisfied. 

We summarize how the time dependent functions £2, £3^ 
A 2 and A3 are determined uniquely and consistently. £ 2 was 
fixed by requesting the normalization condition to be consis- 
tent with the gauge invariance of the initial state, while £3 was 
fixed from the IR regularity of the two point function. As 
presented in Eq. d57| i, to ensure the consistent commutation 
relations, A 2 and A3 are also fixed once £ 2 and £ 3 are given. 
In this paper, we have not derived the gauge-invariance con- 
dition d86l l but we just postulated it. The above discussions, 
however, have proven that this condition can be imposed con- 
sistently by choosing the homogeneous solution appropriately 
inCn- 



10 



V. CONCLUSION 

We presented, in the standard single field inflation model, 
one example of the calculation of a genuine gauge-invariant 
quantity, i.e., the two-point function of 9 R, which is the spa- 
tial curvature perturbation on a </)=constant hypersurface with 
its arguments specified in terms of the geodesic normal coor- 
dinates. We showed that, taking an appropriate initial vacuum, 
the two-point function for 9 R no longer yields IR divergences 
at one-loop order. It would be also possible to extend our argu- 
ment to higher orders in loops and also to the general n-point 
functions. The quantities that are compared with actual obser- 
vations like the fluctuation in the Cosmic Microwave Back- 
ground should also be such genuine gauge-invariant quanti- 
ties. Hence, our result strongly indicates that such quantities 
are also IR regular for the standard single field inflation model. 

In the global gauge that we used in this paper the resid- 
ual gauge degrees of freedom were not fixed. The residual 
gauge degrees of freedom include the overall spatial scale 
transformation corresponding to a constant shift of £ in the 
6<f> = gauge, which is the origin of the IR divergences. To 
remove IR divergences, hence, we had to impose the invari- 
ance of quantum states in the direction of the residual gauge, 
which requests additional gauge invariance conditions, such 
as Eqs. d79l i and d82l i. on the choice of the initial quantum 
state. The condition d79l l restricts the mode function for the 
interaction picture field and the condition ( T82b restricts the re- 
lation between the Heisenberg field and the interaction picture 
field. In the present paper, we derived these conditions, requir- 
ing, instead of the gauge-invariance itself, that the possibly di- 
vergent terms should be written in the form of total derivatives 
under the momentum integral. While the adopted iterative so- 
lution to the Heisenberg equation, that is used to obtain these 
gauge invariance conditions, may not be general enough. In 
this sense, it may be possible to find other vacua that are reg- 
ular against IR contributions. Namely, the gauge-invariance 
condition (=regularity condition) may not uniquely determine 
the vacuum state in the inflationary universe. Therefore what 
we have derived is not a necessary condition but a sufficient 
condition for the IR regularity. We also expect one can spec- 
ify the initial quantum state directly from requirement of the 
gauge-invariance, but this point has not been clarified yet at 
all. We leave these issues for future work. At the leading order 
in the slow-roll approximation, the condition d79l ) is automat- 
ically satisfied in the scale-invariant (Bunch-Davies) vacuum. 
However, at the higher order in the slow-roll approximation, 
this condition looks quite non-trivial. Our preliminary anal- 
ysis tells that this vacuum state seems to coincide with the 
one naturally obtained by using the usual ie prescription. We 
would like to report on this point in our forthcoming publica- 
tion JH]. 

In contrast to the case of the global gauge, if we fix all the 
residual gauge degrees of freedom and quantize only physical 



degrees of freedom, we need not to restrict the initial quan- 
tum state by imposing additional gauge invariance conditions. 
This would provide another way of quantization that also 
yields no artificial divergences. In our previous work I2H1 . 
following this direction, we tried to fix the residual gauge de- 
grees of freedom, including the overall scaling of spatial co- 
ordinates, by imposing appropriate conditions. (We refer to 
this gauge as the local gauge 12111 .) If we perform the canon- 
ical quantization in the local gauge and give the initial state 
there, the local gauge conditions ensure the regularity of the 
IR corrections, without restricting the initial quantum state. It 
is, however, not so trivial to perform canonical quantization 
in the local gauge, because of additional conditions to fix the 
residual gauge degrees of freedom. Furthermore, if we set ini- 
tial vacuum in the local gauge, it would be in general breaks 
the invariance under spatial translation. In our previous work 
I2H1 . We therefore chose the initial vacuum state in the global 
gauge and then linearly transformed the interaction picture 
field from the global gauge to the local one. The truth is that 
in this gauge transformation there appear the non-linear con- 
tributions that can cause IR divergence. These contributions 
were not taken into account in Ref. 1I2TI1 and the expression in 
the local gauge was still including the divergent terms. One 
can say that this is due to the lack of the gauge invariance in 
choosing the initial state, because the IR divergence actually 
disappears at the leading order in the slow-roll approximation 
if we take the Bunch-Davies vacuum, which is invariant under 
the scale transformation. To guarantee the IR regularity also 
at higher orders in this approximation, we need to adapt the 
gauge-invariance condition as is done in this paper 3 . 

In this paper we have solved the Heisenberg equation to the 
second order in the slow roll approximation. Up to this order, 
we showed the presence of a gauge invariant initial quantum 
state that is free from IR divergences. But it looks quite non- 
trivial to extend our results to the higher order in the slow roll 
approximation. For the complete understanding of IR issue in 
the standard single field inflation model, it is also definitely 
necessary to prove the existence of such an initial quantum 
state without relying on the slow roll approximation. 
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Appendix A: Third-order gauge transformation 

In this appendix, we study the change of variables between the flat gauge with S R = and the comoving gauge with 8<f> = 0. 
We denote the gauge transformation from the flat gauge to the comoving gauge as = e £ «x M where £ p is the vector field 
£ p = (a, f3 % ). In the present paper it is not strictly necessary to transform the spatial coordinates, but we do this for future 
application. For the same reason, we also leave the terms that do not yield IR divergences here. 

In accordance with Ref. l44ll . we use the cosmic time t and denote its derivative by a dot. From the condition that this 
transformation makes the scalar field perturbation to vanish, we obtain 

= i P + £^ ( f> + ^ + ^Ll((b + ^ + ^Ll^ + ip) + --- , (Al) 
It is enough to calculate the perturbed expansion up to the third order. Using Eq. ( lAlb . the time shift a is solved at each order as 

f>ai = - T p =: Cn , (A2) 
pa 2 = JLtpV = \(d P + £ 2 )C , (A3) 

P«3 = ^Cnd 2 p Cn + \Cn {dpCnf + ^CudpC + ^(2e 2 + e 3 )(* + ^/Wn , (A4) 

where, (3i = is presumed. In both flat and comoving gauges, the traceless part of the spatial metric is requested to satisfy 
the transverse conditions. To maintain these transverse conditions, we also need to change the spatial coordinates at second 
order. The spatial component of the metric then transforms as 

e 2C [e S % = [e S % + (<* + a 2 + a 3 )e^d t (e 2 " [e%) + P k d k [e 5 % + 2 {d (i $ + [e% )k 

+ \e- 2p { ai + a 2 )d t {(«! + a 2 )d t (e 2p [e^] .J } + V 2 "^ {a x d t {axpe^)} <% + ?%, (A5) 

where Jfy is defined as 

Jfy := 2e" 2p {9 ( ,ai + d {l a 2 ) N n + aie~ 2p 5 t {d^N^ + e 2p <%% 2 } + ^d k pa x 8 i:j + 2 ( 6a 1 <%% 2 
+e~ 2p d {l a 1 |-9j)Q!i - 2<9 i) a 2 - 2d ]) a 1 SN + aiNj) + diiVj) + e^fi^A 
—aiduaidj^ai — ceidiCtidjai . (A6) 

Taking the trace part of Eq. JA5b . we obtain 

C = p{a 1 +a 2 + a 3 ) + {p{a\ + a 2 )} 2 + ^(ai + a 2 )d t {p{ai + a 2 )} + ^e~ 2p a 1 d t {ai<9 t (pa;ie 2p )} 
-(pa,) 2 + i(M) 3 - \P {[e 6 % - [e s %] + \*id t 5«[e s % + l -K 

+ ^(#4 + d^ 3 3 ) [e*] i . - pax ^2pa 2 + ^A, 2 + a ie ~ 2p d t (pa ie 2p ) + ^ , (Al) 

where we defined J{ := S^Oiij. For our purpose, it is sufficient to consider the gravitational wave perturbation up to the second 
order. Neglecting the third-order terms, the transformation of the transverse traceless tensor is given by 

Hi = *7« + Ud p 8% + (Si k S/ + S, k 5 t l - |m W ) (dkU^Ni + d k fr, 2 - lt^ dk C n diCr^ ■ (A8) 

Now, it is clear that Sjij agrees with Sjij, after we neglect the terms that are irrelevant to IR divergences. Using Eqs. ( IA2l i- (IA4l i 
and Eq. (IA8b . Eq. (IA71 > is rewritten as 

C = Cn + \d p £ + \e 2 £ + ±£d 2 p ( n + le 2 Cnd P £ + -^e 2 (e 2 + 2e 3 )£ + Ud P Cnf + \mCn - \< n #fo 

+\{di& + diPl) - IsfidtCn^Nj + l^SfidiCnd^ + Ul - 2( n )H ■ (A9) 
3 3 p z 6 p l 6 

Multiplying the spatial derivative d l on Eq. ( IA5I ). /3j )2 is given as a solution of the Poisson equation: 



d 2 A,2 = -^CndpSjij - &>1Ui,2 + \di-K 2 + -^i 



d k d l ( n d p 6j kl + ( d k d l - -6 kl 8 2 ) Jf fe; , 2 



(A10) 
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where we used <9 ? /3i,2, given by operating d % d^ on Eq. ( IA5b . At the third order, ^3 is obtained in a similar manner. It is notable 
that fi in (n = 2, 3, • • • ) is multiplied by at least one wavenumber vector in the momentum representation. Neglecting the terms 
that are irrelevant to the IR divergences, the curvature perturbation £ is related to as 

C = Cn + \d P C + + \Cd % P Cn + \e*Ud p £ + ^e 2 {e 2 + 2e 3 )£ + ■■■ . (All) 
Here, we also neglected the cubic-order terms with only one gravi ton field Sjjj, since its contribution vanishes in ( 9 R 9 R)±. 
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